Energy as witness of multipartite entanglement in spin clusters 
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We develop a general approach for deriving the energy minima of biseparable states in chains 
of arbitrary spins s, and report numerical results for spin values s < 5/2 (with N < 8). The 
minima provide a set of threshold values for exchange energy, that allow to detect different degrees 
of multipartite entanglement in one-dimensional spin systems. We finally demonstrate that the 
Heisenberg exchange Hamiltonian of N spins has a nondegenerate iV-partite entangled ground state, 
and can thus witness such correlations in all finite spin chains. 
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Entanglement is one of the most striking peculiarities 
of quantum systems and promises to play a crucial role 
in emerging quantum technologies [l|. This has fueled 
the development of theoretical and experimental means 
for its detection in diverse physical systems One of 
the most convenient such tools is represented by entan- 
glement witnesses. These are observables whose expecta- 
tion value can exceed given bounds only in the presence 
of specific forms of entanglement. Macroscopic observ- 
ables such as magnetic susceptibility 0-11] and internal 
energy 041] allow for example to discriminate between 
fully separable and entangled spin states. In qubit sys- 
tems, further inequalities for energy have been derived, 
whose violation implies multipartite entanglement jjjllioj. 
Along the same lines, the measurement of collective ob- 
servables 11] allows to detect multipartite entanglement 
in the vicinity of prototypical quantum states through 
the spin-squeezing inequalities [l2l [l3j . The connection 
of these studies with quantum-information processing has 
however focused most of the attention on entanglement 
between qubits. Limited attention has instead been de- 
voted to multipartite entanglement between composite 
systems of s > 1/2 (pseudo) spins [l4j |. 

In the present paper, we address the problem of de- 
tecting multipartite entanglement in clusters of arbitrary 
spins s [H, [l6| using exchange energy as a witness. We 
develop a general approach for deriving the energy min- 
ima of biseparable states \ipA) <£> IVb) m chains of 
N spins, that exploits the rotational symmetry of the 
system Hamiltonian. This allows to reduce the mini- 
mization problem to calculating the ground states of ef- 
fective spin Hamiltonians within each subsystem A and 
B. Minima derived for A:— spin chains provide in turn 
a set of threshold values for energy, corresponding to k- 
partite entanglement in chains of rik(k — 1) + 1 or rings of 
Uk(k— 1) spins. Analytical expressions of the minima are 
derived for the simplest cases, while numerical solutions 
are provided for s < 5/2, that correspond to prototypical 
models of molecular nanomagnets [17h20| . As a general 
result, we finally demonstrate that the ground state of an 
A-spin chain with Heisenberg Hamiltonian is ./V-partite 
entangled. This implies an energy gap between bisep- 
arable and A-spin entangled states, and the possibility 
of detecting the latter ones by exchange energy, in finite 
spin chains with arbitrary N and s. 

Tripartite entanglement — Tripartite entangled states 



are detected by a three-spin Hamiltonian H123 if their 
energy exceeds the lower bound that applies to bisepa- 
rable states [9(. Here, we seek such bound for H123 — 
Si • s 2 + s 2 • s 3 = H12 + H23, and for a generic biseparable 
state \ipt) ® 1^23): 



Eb s = min {(VilsxlVi) 

|V1)>|V2S) 



(V , 23|s 2 |'023>+ (fe 1 82-831^23) }• 



If we identify the direction of {'023 1^2 1^23) with the z 
axis, the first term in Eq. [T] simplifies to: (H12) — 
(ipi\si : z\ipi}(ip23\s2,z\4>23) , where (fel^^l'fe) > by 
definition. For any given IV23)) the state of s\ that min- 
imizes (H123} is thus given by \mi=— s\), and the prob- 
lem of deriving E bs reduces to finding the state IV23) that 
minimizes: 

(023 1^23 IV'23) = (^23 I - SiS 2: z + S 2 ' S 3 1023), (2) 

i.e. the ground state of the two-spin Hamiltonian i?23- 
In order to derive the energy minima, it is convenient to 
expand IV'23) m the form: 



IV'23) 



M= 




w ; Af\S,M) = 



S 2+S 3 



M=-s 2 - 



S3 



(3) 



where S = s 2 + S3 and M is its projection along z. Each 
real coefficient Pm gives the probability that S has a 
z— projection M (Em ^ m = -0- The normalization con- 

l e %a M reads: 
Given that both the 



dition for the complex coefficients A§ = a s 



M 



EsK 1 ) 2 = 1 (with a% 



operators s 2;Z and s 2 -s 3 commute with S z , the energy ex- 



pectation value can be written as (H23) 
where, 
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The energy expectation value is thus given by an average, 



3 M 

with probabilities Pm , of functions E b ' s that depend on 
disjoint groups of variables A M , each corresponding to a 
given M . This allows to minimize the terms E^ inde- 
pendently from one another, and to identify the overall 



minimum with the lowest E\ 
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-2.481 
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- 3.0 


3/2 


0.749 


0.631 


0.198 


0.0269 






TV 
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-1.933 
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0.676 
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0.0696 


0.00819 




TV 


-8.849 
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-10.0 
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0.120 


0.0232 


0.00244 


7T 


-13.54 


-5.768 


-7.771 


-15.0 



TABLE I: Minima £f s of //123 for biseparable states |^>i) (g> \1p2z) and corresponding coefficients as and as- The ground state 
energies Eq of -H123 are also reported, as well as E12 = — s/(a, a) and £23 = <?(a, a). 



r 3,M 



on the variables Ag 1 is derived 



The dependence of E' b 

as follows. The first contribution in Eq. [4]is proportional 

to: f u = (s 2 , z ) = T.s,A A s)*( A ¥){S,M\s 2 , z \S',M). 
Here, the matrix element can be expressed in terms of 
the Clebsch-Gordan coefficients (S, M\s 2<z \S', M) = 
£ ma (S, M\m 2 ,m 3 )(m 2 ,m 3 \S' , M)m 2 (with m 3 = M — 
m 2 ). The second contribution in Eq. [4] is instead diag- 
onal in the basis \S,M), and reads: gu = {s 2 ■ S3) = 
Es( a ¥) 2 [S(S + 1) - s 2 (s 2 + 1) - s 3 (s 3 + l)]/2. 

In order to analytically minimize - for s < 3/2 - the 



function E^ subject to the normalization constraints, 
we apply the method of Lagrange multipliers. The sta- 
tionary points of the Lagrange function A M (Ag 1 , A) = 
— 1], are identified by the equa- 
dk M /da^ = dA M /d\ = 0, for 
In all the cases considered below, 



Et M + \[Es(4 4 ) 2 



dA M /da¥ 
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tions v±\.mi uUj s 
< S < s 2 - 



S3- 
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the lowest minima correspond to M 
We shall thus refer only to this subspace, and omit the 
apices M from the notation. Besides, we focus on the 
case of identical spins. 

In the s = 1/2 lower bound for (H123) in the 

absence of tripartite entanglement has already been de- 
rived by different means [9(. Here we show that such 
value actually corresponds to a minimum, and derive the 
corresponding biseparable state. The dependence of Ef s 
on the parameters as and as is given by (see Eq. 0J: 
/m = — a oai cos(ao — a\) and gm = (— 3ag + af)/4. As 
far as the phases as are concerned, E% s is minimized by 
«i — «o = 7T. The remaining conditions give rise to the 
energy minimum Ef s = — (l + V5)/4, that coincides with 
the lower bound derived in Ref. [9j. The corresponding 
biseparable state is given by: 



ing biseparable states (Table fl]), that have been obtained 
through a conjugate gradient algorithm (22[. 

The comparison between the different spin values 
shows that the relative weight of the singlet state (ao) 
decreases with increasing s, as well as the ratio between 
the energies of the entangled and unentangled spin pairs 
(E 23 /Ei 2 ). In all cases, the inequality (i/123) < E% s im- 
plies tripartite entanglement in the three-spin system. 
The criterion becomes (H) < n 3 E^ s for any H that can 
be written as the sum of n 3 three-spin Hamiltonians, such 
as chains of 2ri3 + 1 spins or rings with 2n3 . 

Quadripartite entanglement — We consider the ex- 
pectation values of the four-spin Hamiltonian -f/1234 = 
Si • s 2 + s 2 ■ s 3 + s 3 ■ s 4 , corresponding to the biseparable 
states 1^22) = ^12) ® |-034) and \ifjf 3 ) = \tpi) (g) ^234) • In 
the former case, we compute: 



£22 = , m i n , {(^lalsi • S 2 |V>12) + (^34183 • s 4 |V>34 

+ (^12\s 2 , z \^l 2 )(^ 3 4\s 3}Z \tp 3 4)} , 



(8) 



where z is defined as the direction of {^ 3 i\s 3jZ \ip 3 4) ■ 
The states \i>\ 2 ) and ^34) are expanded in the bases 
\S = S 12 ,M = Mia) and \S' = S 43 ,M' = M 43 ), re- 
spectively. For 1^12), we use the expression in Eq. [3] 
(and replace the indices 23 with 12). Similarly, ^34) 
is expressed as: |^ 34 ) = Ea/' VQm' Es B ¥' I s "' M ')' 
with Bff = bg'eVS', Em>Qm> = Es'H-) 2 = 1- 



Being both M and M' good quantum numbers, E 22 = 
Em Em' PmQm'E^ 2 MM , where: 



Et 2 MM = 9M(A M )+f M (A M )f M ,(B M ' 



-g M >(B M ) (9) 



d = 1/2 + 1 



1/2 



1 = (l/2 - 1/V5 



1/2 



(g) and the functions /m and om coincide with those re- 



We proceed in the same way in the case s = 1, where 
the expression of energy is given by: Jm = 2ai(aov2 + 
a 2 )/V3 and gu = — 2og — a\ + a 2 . Here, the conditions 
cts+i—as = 7T, derived from d Km /das = 0, have already 
been included. The analytic expression of the energy 
minimum is: 



E bs 



5/2 cos(y/3) + V3sin(</V3) } , (7) 



where tp = arccos[l/(10vl0)]- 

For the spin values s = 3/2, s = 2, and s — 5/2, we 
directly report the energy minima, and the correspond- 



ported in the previous section. The energy E^ 2 MM 
is minimized numerically by the conjugate gradient ap- 
proach as a function of a A/ and b M , while the minimiza- 
tion with respect to a M and (3 M is straightforward. The 
minimum EL is then identified with the lowest E^ M : 



p4 _ „• p4,MM',-M 

ii/99 = mm fion a 
22 M,W 



a M ,b M '. 



3M' 



(10) 



For all values of s, the lowest minima belong to the sub- 
space M = M' = 0. The minimum of Ef 3 is instead 
identified with the ground state energy of the three-spin 
Hamiltonian H 23 a — ~s\s 2 _ z + s 2 ■ s 3 + s 3 ■ S4, which 
belongs, in all the considered cases, to the subspace 
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- 


- 1.500 


- 1.190 


-1.780 


-2.366 


-2.697 


-3.244 
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0.921 0.387 0.0418 






- 4.051 


- 3.828 


-5.343 


-6.771 


-8.133 


-9.537 


3/2 


0.775 0.607 0.171 


0.0281 




- 8.131 


- 7.957 


-10.90 


-13.75 


-16.56 


-19.39 


2 


0.687 0.669 0.278 


0.0602 


0.00649 


-13.74 


-13.59 


-18.46 


-23.24 


-27.99 


-32.75 


5/2 


0.627 0.669 0.359 


0.134 


0.110 < 10 -4 


-21.18 


-20.71 


-28.03 


-35.23 


-42.42 


-49.62 



TABLE II: Left: Minima E\ 2 and Ef 3 for biseparable four -spin states V12) ® V34) and |V>i) ® IV234), respectively. The 
states corresponding to the former partition are given by the displayed values of as, and by: bg = as, as+i — as = and 
Ps'+l — fis 1 = 0. Right: Energy minima Ej^ of iV-spin systems. 



with M = s. The energy minima and the correspond- 
ing states are reported in the left part of Table HU We 
note that the bipartition IV'12) ® IV34) always gives lower 
minima with respect to IV13) = IV'i) ® ^234): there- 
fore, E^ s — E\ 2 . For the four-qubit system, the ex- 
pectation value of energy is minimized by the dimer- 
ized state @. This is not the case for s > 1/2, where 
the coupling between S2 and S3 induces a significant ad- 
mixture with states of higher S and S'. Besides, the 
energy is minimized by the state with (s2,z) — —(53,2) 



(a = b° and $° s , +1 - 0°, 



"s+i 



7r). We thus 



conclude that, for all the considered spin values, the in- 
equality (-H1234) < Ei s implies quadripartite entangle- 
ment in the four-spin system. The criterion generalizes 
to (H) < n^E^ s for any H that can be written as the sum 
of n4 four-spin Hamiltonians, such as chains of 3«4 + 1 
spins or rings with 3ri4. 

N-partite entanglement — For larger spin numbers N, 
the analytic derivation of the functions Jm and gu be- 
comes cumbersome, and a fully numerical approach is 
preferable. Given a partition of the spin chain in two sub- 
systems, A and B, consisting of Na and Nb — N — Na 
consecutive spins, the Hamiltonian can be written as 
H = Ha + H B + H A b, where H A = E^i _1 s * ' s *+i> 



--N A + 1 



Si ■ Sj+i, and H A b = $n a • s Na+1 . The 

= |Vu) ® IV'b) 



energy minima for biseparable states 
are: 



E 



N 

N A Ni 



= min {(ipA\H A \ipA 

\lpA),\lpB) 



(Vb|Hb|Vs) 

(iPa\sn a \iPa) ■ (*Pb\sn a +i\iPb)} , (11) 



We identify the z direction with that of (^ja\^n a \^a) and 
define: za= (^a\sna,z\^a) > 0, z B = (iPb\sn a +i,z\iPb)- 
Besides, the state IV'b) that minimizes E, 



N A ,N E 



neces- 



sarily has an expectation value (sn a +i) antiparallcl to 
z (and thus z B < 0): any rotation of the subsystem B 
with respect to such orientation would in fact increase 
(Hab), while leaving unaffected (Ha + Hb)- The mini- 
mization can now be split into two correlated eigenvalue 
problems, that consist in finding the ground states of 
H a {zb)=Ha + z b sn a , z and H B (z A ) = H B + z a s Na +i,z- 
The self-consistent solution of the minimization prob- 
lem Eq. QT] is thus represented by the state |V>) = 
IV'a(^b)) ® |V>bM) with: z A = {^ a (zb)\sn a ,z\^a(zb)) 
and z B = (iP b (za)\s Na+ i iZ \4> b (za)}, where \ip A (z B )) 
is the ground state of Ha(zb) and \4>%(za)) that of 



Hb(za)- The corresponding value of energy is given by 
E% aNb = E A (z B ) + E°b(z a ) - z a zb, (12) 

where the last term avoids the double counting of the 
contribution from s^ A ■ sn a +i- The values of the overall 
minima for biseparable states, given by 



E» = min E»* N * 

N A ,N B iV 



(13) 



are reported in the right part of Table |TT] for N < 8. 
For all the considered values of s and N, the partition 
with lowest energy minimum is that with Na — 2. We 
note that for even Na and N B , the qubits only present 
a solution with (Hab) = 0; for s > 1/2, instead, the 
minimum corresponds to the additional solution, with 
finite (Hab)- As in the cases of tri- and quadri-partite 
entanglement, these minima provide a criterion, namely 
(H) < nkE^ s , for the detection of fc-partite entanglement 
in chains and rings with nk(k— 1) + 1 and rif.(k— 1) spins, 
respectively. 

We finally demonstrate the presence of TV-partite en- 
tanglement in the ground state of all spin chains with 
even N. 

Theorem. — The ground state IV'o) of the spin Hamil- 
tonian H = X^fcT 1 s * ' s i+ii with even TV, cannot be writ- 
ten in any biseparable form \iPab) = |Vu) ® IV'b)) and is 
thus iV-partite entangled. 

Proof. — According to Marshall's theorems [23|, |Vo) 
is a nondegenerate singlet state. 

A biseparable state |Vm,b) can only be a singlet 
if S x = (x = A,B). In fact, one can write 
) as a linear superposition of eigenstates of S 2 : 



\1>x 

IV'x) = E s - 

(S 2 ) > Es A ,s 

^2s A ,Sb 



The following inequality applies: 



clcL\ 2 [(S a - S B ) 2 + S A + S b ] > 



r<A r<B 

^Sa^Sb 



'S A ^S 

\ 2 (Sa + Sb), where we make use of: 
^aI^sJ-^sJ^V > SaSb- Therefore, (S 2 ) 
can only vanish C£ A = 5 Sa , and C§ B = S Sb ,q. 

We now prove that the state |Vu) <£> IV'b), with S A /b = 
0, cannot be the ground state of H, by showing that 
-H'IV'ab) has a component |Vab) which is orthogonal to 
|Vub)- To this aim, we write: H A b = s z ,n a s z ,n a +i + 
(s +y N A S-,N A +i +h.c.)/2. We first show that |Vis) = 
Sz.tVaSz.a'a+iIV'ab) is finite and belongs to the subspace 
5*4 1 ft = 1 and Ma/b = 0. In the partial spin sum basis 
[2l|, the state of A reads: (Vu) = J2 a D a \a, Sa, Ma), 
where a denotes the quantum numbers S±, . . . , Sn a -i 



4 



corresponding to the partial spin sums Sfe = Ym=i s fc> 
and Sa = implies Sn a -i — s. The operator s z Na com- 
mutes with all S\ with k < Na — 1- The matrix elements 
of the Na— th spin can thus be reduced to those between 
the states of two spins s: (a', S' A , M' a \s z Na |a, Sa, Ma) = 
$u,a'{S' A ,M l A \s z NA \Q 1 ()). The latter matrix element is 
only finite, and equals — i] s , for S' A = 1 and M' A = 0; 
therefore, s^JVu) = ~Vs E« D a K !> °)> with % = 

[(EL- s ™ 2 )/( 2s + !)] 1/2 > °- The same proce- 
dure can be applied to B, resulting in: s z Na+1 \^ b ) = 

-VsEp Df\P, 1,0). Here |V>b> - T,pDf\^S B = 
0, Ms = 0), and /3 denotes the quantum numbers 
Si,..., S Nb -i corresponding to S k = J2i=i s N+i-i- As 
a result, \^ B ) = rj 2 s Y, a ,p K Df\a, 1, 0) ® |/3,1,0> has 
finite norm, belongs to the subspace Sa/b = 1 an d 
M A /b = 0, and is thus orthogonal to |Vub)- 

We now show that IV'ab) coincides with the compo- 
nent of H\iPab) with Sa/b = 1 an d Ma/b = 0. In 
fact, (ii^4 + Hb)\4>ab) belongs to the SU/b = sub- 
space, being [if x ,S x ,] = for x,x' = ^>-B- The 
states s^s^+iIV'ab) belong instead to the subspaces 
Ma = —Mb = ±1. Therefore, iZ^As) has a finite com- 
ponent IV'ab); an d cannot be an eigenstate of H. 

We finally consider the case where the spins of the 
subsystems are not consecutive. In the simplest case, 
the spins of A are split into two sequences of and 
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Na 2 consecutive spins, separated by the Nb spins of B. 
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In conclusion, we have developed a simple approach for 
deriving the energy minima of biseparable states in chains 
of arbitrary spins s. These minima can be used for detect- 
ing fc-partite entanglement in chains with n^(k — 1) + 1 
and rings with rik(n — 1) spins, respectively. This ap- 
proach has been here applied to spin chains of up to 8 
spins s < 5/2. Finally, we have demonstrated on general 
grounds that the Heisenberg interaction induces N par- 
tite entanglement in the nondegenerate ground state of 
even-numbered chains with arbitrary s. Such entangle- 
ment can thus always be detected by using energy as a 
witness. 
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